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Abstract. The cosmic microwave background anisotropies produced by active seeds, such 
as topological defects, have been computed recently for a variety of models by a number of 
authors. In this paper we show how the generic features of the anisotropies caused by active, 
incoherent, seeds (that is the absence of acoustic peaks at small scales) can be obtained semi- 
analytically, without entering into the model dependent details of their formation, structure 
and evolution. 
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1 Introduction 

Two cosmological scenarios are currently in competition to explain the large scale 
(6 > 7°), fairly flat, spectrum of cosmic microwave background (CMB) anisotropies 
which has been observed by the COBE satellite jjj : inflation (see, e.g. 0j), and active 
seed models (see, e.g. ||). 

Observations of the CMB anisotropy spectrum on smaller scales have been per- 
formed by balloon and ground experiments B, which tend to indicate the presence 
of a peak at 9 ~ 1°. This peak in the spectrum gives at present the lead to the 
inflationary scenario. 

One should however be careful not to exclude too hastily active seeds at this stage. 
Indeed, to start with, the predictions from specific topological defect models are not 
yet very robust, the reason being that they act as a continuous, non-linear source of 
inhomogcneities, and hence are difficult to model. Second, there exist simple models 
of active seeds which "mimic" inflation and reproduce a peak in the spectrum on the 
degree scale |J. Third, including some microphysics in the evolution of the defect 
network modifies the "standard" picture 0. Even if such modellings are arguably 
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unrealistic, they show that the absence of a peak is not the seal of all types of active 
seeds. 

It is also important to be able to predict which active seeds produce secondary 
peaks, like inflationary models, and which do not. As we shall see, the absence of 
secondary peaks is the generic signature of "incoherent" seeds described within the 
"stiff" approximation. More then that the presence of a sharp peak on the degree 
scale, the absence of secondary peaks, which will be probed by the MAP and Planck 
satellite missions [M || (and possibly also by the Boomerang and Archeops |], |l0) 
balloon experiments), will toll the bell of active seeds as the main contributor to CMB 
anisotropies. 

It is clear that the (semi-)analytic calculations of the seed stress-energy tensor 
correlators pioneered by Durrer and collaborators p"l| ] have the advantage, over heavy 
numerical calculations within specific topological defect models, to test the influence 
on the CMB anisotropies of each component of the seed stress-energy tensor. This 
semi-analytic approach has already allowed to grasp some generic features of the CMB 
anisotropies seeded by active sources (see e.g. [|[ [| [ll], |l2[ |ll| 0). 

The object of this paper is to study, within this semi-analytic approach, incoherent 
active seeds and show that, within the stiff approximation, they all lead to CMB 
spectra exhibiting no secondary peaks. 

The paper is organised as follows : in we first recall the some statistical prop- 
erties of the seed stress-energy tensor correlators at large wavelengths [see eq. (^6|)j, 
and write some up-to-now unmentionned relations [see eqns. (|27j)l. We then give the 
new result that any active source described by a test scalar field must obey these 
properties [see eqns. (^-^)]. In §0, we first recall the definition and properties [see 
eq. (|46|)1 of coherent sources (e.g. jl3|, 0), and then explicitely construct for the first 
time a generic (=incoherent) active source as a sum of coherent ones. Finally, in §^|, 
we exhibit the generic behaviour of the CMB anisotropies seeded by active sources. 



2 The statistical properties of incoherent sources 

2.1 The 2-point correlators 

The stress-energy tensor of active sources is a small perturbation added to the other 
cosmic fluids (this is the so-called "stiff approximation" (see e.g. We decompose 

its components 9 M „(a; 1 into their scalar, vector and tensor parts (SVT) as (/i, v = 
0, 1,2,3, rj is the conformal time, x l are 3 cartesian coordinates; space is assumed to 
be flat for simplicity) : 

e 00 = p s , (i) 

e 0l = -(vt + d lV s ), (2) 
Qij = <\J- ■ (i'.j ^,,a) ii- ■ n,\r ■ rv.n; • lC (3) 

Barred spatial vectors are divergenceless; barred spatial tensors are traceless and di- 
ver genceless. 
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We work in Fourier space, the Fourier transform of any function f[x % ,rj) being 
defined as : 

f(k\v) = J e- ik ^f(x\r,)d s x f(x\ V ) = ^ J e ik ^ f(k\ V )d 3 k. (4) 
From ([[III); we therefore have : 

kk 3 -h t3 kAfl s = (\P 13 P U +L k L l \® kh (5) 



2 *(* n i) 



(P^ l l)e kl , (6) 



n y = [p k P l j --p ij p kl )e kl , (?) 



and similar expressions for the other variables, with 



k- k 

Lij = Sij Pij — — kikj . (8) 



The ten components of Q^ v {x l of the active source are ten statistically spa- 
tially homogeneous and isotropic random fields. The statistical properties of those ten 
random fields are described by their unequal time two-point correlators 

(&^(x l ,T])e p a(x'\r]')) =C fl „ pcT {r\Ti,ri'), (9) 

where (...) means an ensemble average on a large number of realisations, and where the 
correlator C^ uprT is a tensor which depends only on r], rj' and r % = x l — x n because of the 
spatial homogeneity of the distribution. The power spectra of the Fourier transforms 
are given by : 

(e;^,77)e pa (fc'W)) = 8(k i -h*)d lu , pa (k i ,ri, rf). (10) 

The spatial isotropy of the distribution now forces the power spectra to be of the form 



Coooo 


= A , 






(11) 


CoOOi 


= ikiBx, 






(12) 


CoOij 


= Co5 i: j + C 2 kikj, 






(13) 


Coioj 


= D 5ij + D 2 kikj, 






(14) 


Coijk 


= i [EikiSjk + Ei(kjSik 


+ kkSij) 


+ E 3 kik 3 k k ] , 


(15) 


Cijki 


= F SijSki + F (S i i ! Sji -\ 


- Sitfjk) ■+ 


- F 2 (k i k :j d k i + k h ki5ij) + 


(16) 




F 2 (kikkSji + hkiSjk - 


f kjkiSik 


+ kjkkSu) + F±kikjkkh, 





where Aq, B\ etc are 14 real functions of rj, rf and the modulus k of the spatial vector 
k\ 

We suppose that the active sources appeared at a definite time so that the distribu- 
tion must be, for causality reasons, completely uncorrelated on scales larger than the 
particle horizon (we assume a standard Big-Bang scenario). Therefore, as stressed e.g. 
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by Turok fll2fl , the unequal time correlators are strictly zero outside the intersection 
of the past light-cones, that is : 

CW(r\77,7/) =0 if r>r 1 + r ) '. (17) 

Property (|l7]) translates in Fourier space into the fact that the equal time power 
spectra are white noise on super-horizon scales (that is for kr\ <C 1). Indeed, because 
the correlators ( |l7j ) have compact supports, their Fourier transforms are C°° in k l . 
Therefore causality forces the fourteen functions Aq, B\ etc to be C°° in k 2 . Moreover, 
since within one horizon volume there are almost no sources, those fourteen functions 
must tend to zero on small scales, that is for kr\ S> 1. When is decomposed into 
its scalar, vector and tensor components according to eqns. (|l|||), then eqns. ( |lC)| ) and 
§§ yield 



V;V s j = D a + k 2 D 2l (18) 
(C% S ) = D P id , (19) 



= -fc(2^i + A: 2 £; 3 ), (20) 

(C^ S ) - -kEiPa, (21) 



(tt* s tt s ) = 3F Q + Ak 2 F 2 + k 4 F 4 , (22) 
(0;) = (Fo + k 2 F 2 )P lJ: (23) 

fl*K.i^ = FoiPikPjt + PuPjk-PijPu), (24) 

and similar expressions for the other correlators (see e.g. [fl4|]), where we have defined 
^p*(k i ,j])p 8 (k H ,j]')) = S(k l — k n ) {p* s p s } etc, and where we have introduced 

tt s = k 2 tl s , ft* = fcfl- , V s = kv s . (25) 

An important property of these correlators is that since (for causality reasons) 
^^2) and (k 4 F4j are generically of higher order in k than Fq, then, for small k [|12[ 

i (Ctt,) ~ \ (f TK) * \ (CC ) ■ (26) 

Similarly, when (k 2 D 2 ) is of higher order in k than Do, and when {k 2 E^) is of higher 
order in k than E\, one has : 

(^4(CC) , (K^)^{vtX)- (27) 

These ratios are due to the geometric properties of the SVT decomposition and 
the non-linear structure of the seed stress-energy tensor. In order to give some insight 
about this phenomenon, we shall present in the next paragraph the example of a scalar 
field. 
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2.2 An example 



We consider a TV-component real scalar field ip A evolving in an arbitrary potential V 
according to the general Klein-Gordon equations : 



dV 



= 0. 



(28) 



This equation describes the evolution of a large class of topological defects. The 
spatial traceless part of the stress-energy tensor (denoted with a subscript ST) does 
not depend on V, and reads (see also flq]) : 



e 



ST 



(29) 



where a sum on the index A is understood, and where * stands for the convolution 
product, so that 



e 



ST 
ij 



= - Xijip(p, n)tp(k - p, n)dp, 



P(i(kj) -Pj)) - -{p.k-p 2 )5 ir 



(30) 
(31) 



We decompose Xy into its scalar vector and tensor parts, denoted respectively by the 
superscripts S, V, T: 



X 



X 



k 



1 



3~ u 



k^ kj 



1 k 3 , 



p 



2(i 



pikj H- kipj 



jjik^ kj 



V T 2 
X ti = P 



~5ij(l - (?) - -kikj(l + (i 2 ) -pipj + p{pikj + ki-pj) 



(32) 
(33) 
(34) 



with the notations 

ki = ki/k , pi=pi/p , (i = pik l . (35) 
Now, the correlators of the spatial traceless part of the stress-energy tensor are : 

^ef/^fc^ef^fc',//)) = j ' x i:i {p,k)x H (p ',£')« (36) 

x (il>*(p,T})ip*(ie - p,rj)ip{jf ,rj')ip(k' -jf,r)') 



Since the source distribution is homogeneous and isotropic, this expression can be 
written under a form similar to eqns. ( ^0[]l^ ) : 



e^*(M)e£T(*',»/)) = s(k-k') (37) 

x / Xij(p,k^Xki(p,k)Q(k,\p-le\,T),Ti')dp. 



A. Riazuelo et al., CMB anisotropies seeded by incoherent sources 



293 



We see that the angular dependance arises only in the tensors Xy. Moreover, in the 
limit k — > 0, one has 

(The scalar, vector and tensor parts are always decoupled.) In addition, we see that 
in the same limit k — > 0, all the correlators ([37|) are proportional to the integral 

/ = J p 6 Q(0,p,7 7 ,7/)dp, (39) 

(which is always non zero, being the average of a quartic term) so that using eqns. (|^-0, 
p7H39| ), one gets : 

n s n s \ = Tj x 3 (n s n s ) =Iax3, (40) 

n l s nA = §x2 (^' 8 \ = /«x2, (4i) 



fc 2 

n «n s 1 = ^x4, (42) 

where a is a numerical factor of order unity. It is easy to show along similar lines that 
the correlators involving the velocities obey eqns. (|27|). 

Several conclusions arise from this simple, although generic, model. The first is 
that, in the long wavelength limit, one finds, as expected, that the scalar, vector and 
tensor correlators of the anisotropic stress ( p2[]2~4| ) are in the ratio 3:2:4. The second 
is that this ratio is only an artefact of the SVT decomposition and of the fact that 
the seed stress-energy tensor is quadratic in i\> : it does not depend on the detailed 
dynamics of ip- However, due to the presence of the k/p terms in the X-j , as well 
as the presence of k in Q [see eqns. (|||HIZ|)], 

the angular dependance for larger k 

will no longer be in the same 3:2:4 ratio : the function Fq defined in eq. (E3) is 

s v 

therefore proportional to /, and the higher order terms in k in X^ and Q generate 
the expressions for Fi and F4. 



3 Active sources as sums of coherent ones 

3.1 Coherent sources 

By definition, a coherent source is such that the correlators (fTl-16) factorise 



C MVp ^(v,i ',*:*) =c^(k\n)c pa (k t ,r)'). (43) 
Such a requirement implies, as shown in details in |l4|] , that : 

p s = VHfi(u, V )e(k) , P 8 = VHf 3 (u,r))e(k), (44) 

V s = VHuf 2 (u,n)e(k) , = VHuf 5 (u, rj)Si(k), (45) 

tt s = Vnu 2 U(u,n)e(k) , W* = VHti 2 / 6 (M,r/)e 4 (fc), (46) 

ft*,- = Vnu 2 f 7 (u,r,)E tJ (k), (47) 
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where TL is the comoving Hubble parameter, where u = k/TL. All the /, are arbitrary 
functions behaving as u° when u — > 0, and decaying to when u — > oo. When the fa do 
not depend on ry, the sources have a scaling behaviour. The e, ej, e~ij are independant 
complex random variables which can be defined in such a way that 

(e(k)e*(k')) = (e\k)e*(k')) = (l« = S(k-k'). (48) 

An immediate consequence of (@-||) is : 

, (CC) , (n*^s)-o(kyH 3 ) (49) 

and that, contrarily to the general case (p6|), the anisotropic stress correlators of 
coherent sources are not in a definite ratio for small k. The question to be asked is 
whether a sum of coherent sources (with for example anisotropic stresses correlators of 
order k A ) can lead to a generic incoherent source (with anisotropic stresses correlators 
of order k°). 



3.2 Coherent decomposition of a generic source 

As already stressed by several authors [jl6], [l?], an incoherent source can formally 
be decomposed into a sum a coherent eigenmodes, that is its correlators can be written 
as : 

<Wr(fr, V, rf) = £ A«cW (fc, r7)c«(fc, ?/). (50) 
W 

By coherent decomposition, we mean that the eigenmodes Cpi, behave at low k like the 



coherent sources of eqns. fl44f|47] ) . In a given topological model, these eigenmodes are 
extracted from the C fivpa (k, rj, rj') given by the numerical simulation of the network |l6| , 

In the semi-analytic approach adopted here, one can conversely construct an active 
source by summing several functions fj (for the tensorial part) behaving as k° but in 
such a way that their sum behaves as k~ A so as both eqns. ( p4j ) and ( p7| ) are satisfied. 
We show here a simple example of how this can be realized. 

Let us consider a set of N functions f n {k) obeying causality and white noise con- 
straints [i.e. fn(k) oc k° when k — > 0, and /„ — > when k — > oo]. For simplicity only, 
we consider 

f n (k) = A n Y(k n - k), (51) 

A n and k n being two sets of positive real numbers. We also impose for simplicity 
that Vn, fc„+i < k n . Let us consider then 

F N = Y,fn- (52) 

a 

It is easy to see that 

ffflKfco) = y, /n(feo) = /]An with n = max{n|fe n > k }. (53) 



n=0 
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We then consider a function F which is decreasing on R + and tends towards oo around 
(for example, F(k) = k~ a ). It is straightforward to see that by choosing 

A n = l/y/N , k n =F- 1 (n/^N), (54) 

the set F/v converges towards F. We have therefore explicitely built a function behav- 
ing as l/k a with an (infinite) sum of functions behaving individually as k°. Coming 
back to the correlators, the cost is that the causality constraints fade away as we in- 
crease the number of functions in the sum since x n — > when N — * oo. However, the 
causality constraints need not be satisfied by the coherent eigenmodes, only the sum 
requires causality to be satisfied. 

The fact that we have taken a Heaviside function instead of a more regular function 
is not important. A set of gaussians of various widths would give the same result at 
the leading order in k. The fact that we use a finite truncated sum is unimportant as 
well, as long as the eigenmodes we neglect influence scales much bigger than today's 
Hubble radius only. 



4 Results and conclusion 

We considered here tensor modes only, and we computed their contribution to the 
CMB anisotropies by means of a Boltzmann code developped by one of us (A.R., 
see [|l9| for the detailed equations). We took a sum of coherent and scaling modes 
such that 

"•(«) m 2 nr. , ,2.2 



/WuV^nexp(-fe>72)ei j! (55) 
with A n and k n given by eqns. (|54|), such that 

(^ij in) ^,(n)) a k °- ( 56 ) 

n 

The results are shown on fig. [l| We see that the sum converges, and that some of the 
modes have roughly the same amplitude in the region of the Doppler peaks, but not 
the same phase : the Doppler peak structure is washed out. 

The disappearance of Doppler peaks comes from the fact that although each coher- 
ent eigenmode gives an oscillatory contribution to the CMB anisotropy, their sum does 

= s '( n ) 

not. Indeed, as n increases, the extension in u of decreases. This implies that 

the metric is perturbed for shorter and shorter duration hence on larger and larger 
scales. 

The fact that one has to consider a sum of several eigenmodes comes from eq. (|37|), 
which is a convolution product, hence is not factorisable as a product of functions 
c(k,r))c(k,r)'). However, it is possible that one eigenmode dominates the sum, in 
which case the Doppler peak structure can partially remain, as seems to be the case 
in the large N model pQ|. 



Finally, the convolution product appearing in eq. (36) comes from the fact that any 
stress-energy tensor component is (a least) quadratic with respect to the perturbation 
(i.e. the field ip), therefore any correlator is at least quartic with respect to the field. 
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Fig. 1 Analytic construction of the tensorial contribution of incoherent active sources to 
the CMB anisotropies. The incoherent source (solid line) is a sum of 12 coherent eigen- 
modes (dashed lines). The correlators of the coherent eigenmodes behave as k , whereas the 
correlator of the incoherent source behaves as A: . Vertical units are arbitrary. 



In the most optimistic case, by using Wick's theorem (this is possible when one studies 
the large N model, see ^0|), this quartic contribution can be reduced to a quadratic 
one, but even in this case, the remaining formula eq. (|3^) is not factorizable as a 
product of two quantities. 

On the contrary, in inflationary scenarios (and, more generally, in any scenario pro- 
ducing an initial power spectrum of fluctuations), the cosmic fluids (baryons, neutrinos, 
etc) are described as coherent sources. Indeed their correlators are only quadratic in 
the perturbations, because there exists a contribution of the fluids to the background. 
For example, one has (forgetting here the metric perturbation) : 

8T 00 oc p(rj)8(k,rj), (57) 
STij oc P(r))TTij(k,ri), (58) 

and so on, and solving the equations of evolution, one will find that the perturbations 
are linear with respect to the initial conditions (as in the case of the large N model, 
but here, the stress energy-tensor is also linear with respect to the perturbations) : 

8(k,T,) = £*(a)(M)*(«)(M*), (59) 
(a) 

where Xt a ) stands for all the 8, V, it of the different species which take part in the 
initial conditions (at r\ — rf). Then if one calculates the correlators of the stress- 
energy tensor components, one obtains a quadratic dependance with respect to the 
perturbations, which can be cancelled by using the correlators of the initial conditions. 
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For example : 

(6T* m (k,r,),6T o(k',v')) 
oc PWp(^)^(a)(fc^)i ;l ( 6 )(fc',r/)(l ( * a) (fc,r / *)l (b) (fc',r ? *)) 

(a), (6) 

oc £ p(v)pW)F {a) (k )V )F {b) (k\r,')S(k-k')C {ab) (k). (60) 

(a), (6) 

Moreover, all the perturbations generally are supposed (but this is not always the case, 
see |2^]) to depend on only one random variable so that the sum in the last expression 
is reduced to one component. We also see that one will not get any convolution product 
of the form (j>, r])F^ (\k — p|,r/)j an d therefore the correlators will factorise. 

In conclusion, topological defects are incoherent because their stress-energy tensor 
is quadratic with respect to the field, whereas for standard fluids (i.e. fluid which have 
a contribution to the background) it is only linear. Therefore, the fact that defects 
are incoherent is ultimately linked with the stiff approximation. 

It is a pleasure to thank Martin Kunz and Jean-Philippe Uzan for helping us to clarify some 
points about the O(N) model, and Ruth Durrer for numerous, fruitful, enjoyable discussions. 
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